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Abstract. We consider the problem of finding minimum-link rectilinear paths in rectilinear polygonal 
domains in the plane. A path or a polygon is rectilinear if all its edges are axis-parallel. Given a set V 
of h pairwise-disjoint rectilinear polygonal obstacles with a total of n vertices in the plane, a minimum- 
link rectilinear path between two points is a rectilinear path that avoids all obstacles with the minimum 
number of edges. In this paper, we present a new algorithm for finding minimum-link rectilinear paths 
among V. After the plane is triangulated, with respect to any source point s, our algorithm builds an 
0(n)-size data structure in 0{n -I- h log h) time, such that given any query point t, the number of edges 
of a minimum-link rectilinear path from s to t can be computed in O(logn) time and the actual path 
can be output in additional time linear in the number of the edges of the path. The previously best 
algorithm computes such a data structure in 0(n log n) time. 


1 Introduction 

Paths with few turns have applications in a variety of areas, including VLSI design, computer vision, 
cartography, geographical information systems, robotics, computer graphics, image processing, and 
solid modeling [2112117122125129130134136] . Finding paths with few turns (or minimum-link paths) 
has received much attention, e.g., |1I2I9I1()I11I14I16I19I23I24I26I27I3()I32I33I34I35| . In this paper, we 
enrich the literature by presenting a new algorithm for finding the minimum-link rectilinear paths 
among rectilinear polygonal domains in the plane. 

Given a set V oi h pairwise-disjoint polygonal obstacles with a total of n vertices in the plane, 
the plane minus the interior of all obstacles is called the free space. The link distance of a path is 
defined to be the number of edges (also called links) in the path. A minimum-link path (or min-link 
path) between two points s and t is a path from s to t in the free space with the minimum link 
distance. The min-link path query problem is to construct a data structure (called link distance 
map) with respect to a given source point s, such that for any query point t, a min-link path from s 
to t can be quickly computed. In the following, we say a link distance map has the standard query 
performance if given any t, the link distance of a min-link s-t path can be computed in O(logre) 
time and the actual path can be output in additional time linear in the link distance of the path. 

A polygon (or path) is rectilinear if all its edges are axis-parallel. P is a rectilinear polygonal 
domain if every obstacle of V is rectilinear (e.g., see Fig. [T]). The rectilinear version of the min- 
link path/query problem is to find min-link rectilinear paths in rectilinear polygonal domains. 

* A preliminary version of this paper will appear in the Proceedings of the 42nd International Colloquium on 
Automata, Languages, and Programming (ICALP 2015). 
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Fig. 1. Illustrating a rectilinear min-link path in a rectilinear polygonal domain. 

Rectilinear polygons are commonly used as approximations to arbitrary simply polygons; and they 
arise naturally in domains dominated by Cartesian coordinates, such as raster graphics, VLSI design 
(as mentioned earlier), robotics, or architecture. 

1.1 Previous Work 

Linear-time algorithms have been given for finding min-link paths in simple polygons [16I18I33I34I35] . 
The link distance map can also be built in linear time [33134135] for simple polygons, with the 
standard query performance. For polygonal domains, the problem becomes much more difficult. 
Mitchell, Rote, and Woeginger [27j gave an 0(n^a(n) log^ n) time algorithm for finding min-link 
paths, where a{n) is the inverse Ackermann function; a link distance map with slightly more con¬ 
struction time is also given in m- As shown by Mitchell, Polishchuk, and Sysikaski [26|, finding 
min-link paths in polygonal domains is 3SUM-hard. 

The rectilinear min-link path problems have also been studied. For simple rectilinear polygons, 
de Berg m presented an algorithm that can build an 0(n)-size link distance map in 0{n log n) time 
and 0{n) space, with the standard query performance. The construction time was later reduced to 
0{n) time by Lingas, Maheshwari, and Sack [23], and by Schuierer [32] . 

For rectilinear polygonal domains, Imai and Asano m presented an 0(n log n) time and space 
algorithm for finding min-link rectilinear paths. Later, Das and Narasimhan |9] described an im¬ 
proved algorithm of 0{n log n) time and 0(n) space; Sato, Sakanaka, and Ohtsuki [3l] gave a similar 
algorithm with the same performance. Recently, Mitchell, Polishchuk, and Sysikaski [26] presented 
a simpler algorithm of 0(n log n) time and 0(n) space. Link distance maps of 0(n)-size can also be 
built in 0(n log n) time and 0(n) space [9126] . As shown in [9124] . the problem of finding min-link 
rectilinear paths in rectilinear polygonal domains has a lower bound of f2{n+hlog h) on the running 
time. Therefore, the algorithms in |!)l2(il.;ij are optimal only if h = 0{n). 

However, since the value h can be substantially smaller than n, it is desirable to have an 
algorithm whose running time is bounded by 0(n -|- f{h)), where f{h) is a function of h. 

In addition, the C-oriented version of the min-link paths problems have also been considered 
|1|18|26] , where the edges of the paths are allowed to have C different directions. Our rectilinear 
version is essentially a special case of the C-oriented version with C = 2. 

1.2 Our Results 

We consider the rectilinear min-link paths in a rectilinear domain V. After the free space of V 
is triangulated, our algorithm builds a link distance map in 0{n + hlogh) time and 0(n) space, 
with the standard query performance. The triangulation can be done in 0(n log n) time or 0{n + 
hlog^^^ h) time for any e > 0 |3]. Hence, our result is an improvement over the previous 0(nlogn) 
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time algorithms |9l26l31j . especially when h is substantially smaller than re, e.g., if /i = 0(re^ for 
any 5 > 0, our algorithm runs in 0(re) time. 

1.3 Algorithm Overview 

Our idea is to combine Das and Narasimhan’s algorithmic scheme [9] and a corridor structure of 
polygonal domains. The corridor structure and its extensions have been used to solve shortest path 
problems, e.g., |bl7l2()l28j : however, to the best of our knowledge, it was not used to tackle min-link 
path problems. The corridor structure partitions the free space of V into 0{h) corridors and 0{h) 
“junction” rectangles that connect all corridors. 

The algorithm in [9] (which we call the DN algorithm) sweeps the free space, from the source 
point s, to build the map. The sweep is controlled in a global way so that the time is bounded 
by O(relogre). This global sweeping on the entire free space restricts the DN algorithm from being 
implemented in 0{n + hlogh) time because each operation takes O(logre) time and there are 0(re) 
operations. Using the corridor structure, our algorithm avoids the global sweeping on the entire free 
space. When the sweep is in junction rectangles, we control the sweep in a global way as in the DN 
algorithm. However, when the sweep enters a corridor, we process the corridor independently and 
“locally” without considering the space outside the corridor. Since a corridor is a simple polygon, 
we are able to design a faster algorithm for processing the sweep in it. When we finish processing 
a corridor, we arrive at a junction rectangle. Next, we pick an unprocessed junction rectangle that 
currently has the smallest link distance to s to “resume” the sweep. This is somewhat similar to 
Dijkstra’s shortest path algorithm. In this way, there are only 0{h) operations that need to be 
performed in logarithmic time each. 

To achieve the 0{n + hlogh) time bound, we have to implement the algorithm in a very careful 
manner. For example, we need an efficient algorithm to compute link distance maps in corridors 
(i.e., simple rectilinear polygons). Although efficient algorithms are available for computing the link 
distance maps in simple polygons, e.g., [191181231,32] . they are not suitable for our main algorithmic 
scheme, which requires an algorithm with some special properties. Specifically, let C be a corridor. 
Suppose there are hi pairwise-disjoint segments on a vertical edge di (called a “door”) of C and 
these segments are considered as “light sources”, stored in a balanced binary search tree. We want 
to build a link distance map in C, and obtain a balanced binary search tree storing the light sources 
(let /i 2 be their number) going out of C through another vertical edge d 2 (another “door”) of C. 
Our goal is to design a “corridor-processing” algorithm for the above problem with running time 
0{m + {hi — h 2 + 1) log hi), where rre is the number of vertices of C. We derive such an algorithm, 
which might be of independent interest. 

The rest of the paper is organized as follows. In Section [21 we define notation and review the DN 
algorithm [9]. In Section [3l we introduce the corridor structure and present the main scheme of our 
algorithm while leaving our algorithms for processing corridors in Section [H Section [5] concludes 
the paper. 

2 Preliminaries 

For simplicity of discussion, let TZ he a large rectangle that contains all obstacles of V and let 
iF denote the free space of V in TZ. Note that our algorithm can handle the case where TZ is an 
arbitrary rectilinear polygon, but for simplicity of discussion, we consider the case where 7^ is a 
rectangle. We assume IF has been triangulated. Let s be a given source point in TF. For ease of 
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exposition, we make a general position assumption that no three vertices of P U {s} have the same 
X- or y-coordinate; the assumption can be lifted without affecting the performance of the algorithm 
asymptotically although more tedious discussions would be needed. In the following, paths always 
refer to rectilinear paths in For any point t in a path from s to t is also referred to as an s-t 
path. An edge of a path is also called a link of the path. 

Consider any point t ^ T. An s-t path tt is called a horizontal-start-vertical-end path (or h-v- 
path for short) if the first link of vr (i.e., the edge incident to s) is horizontal and the last link of (i.e., 
the edge incident to t) is vertical. The h-h-paths, v-h-paths, and v-v-paths are defined analogously. 
To make it consistent, if vr is an h-h-path of k links, we also consider it to be an h-v-path of /c + 1 
links (i.e., we enforce an additional edge of zero length at the end of the path), and similarly, it is 
also considered to be an v-h-path of A: + 1 links and a v-v-path of /c + 2 links. A min-link h-v-path 
from s to t is an h-v-path from s to t with the minimum number of links. The min-link h-h-paths, 
min-link v-h-paths, and min-link v-v-paths are defined analogously. To find a min-link s-t path, our 
algorithm will find the following four s-t paths: a min-link h-v-path, a min-link h-h-path, a min-link 
v-v-path, and a min-link v-h-path. Clearly, among the above four paths, the one with the minimum 
number of links is a min-link s-t path. 

To answer the min-link query, our algorithm will compute four link distance maps of 0(n) 
size each: an h-h-map, an h-v-map, a v-h-map, and a v-v-map, defined as follows. The h-h-map is 
a decomposition of J- into regions such that for any region R, the link distances of the min-link 
h-h-paths from s to all points in R are the same. The other three maps are defined analogously. In 
addition, we build linear-size point location data structures on these maps in linear time. 

With the above four maps, for any query point t, we determine the region containing t in each map 
and the one with the smallest link distance gives our sought min-link s-t path distance. 

The vertical visibility decomposition of F, denoted by VD{F), is obtained by extending each 
vertical edge of the obstacles in V until it hits either another obstacle or the boundary of TZ (e.g., 
see Fig. [2]). We call the above edge extensions the diagonals. We consider the source s as a special 
obstacle and extend a vertical diagonal through s. Since F has been triangulated, VD{F) can be 
obtained in 0(n) time |3|4|5|15] . In VD{F), F is decomposed into rectangles, also called cells. Due 
to our general position assumption, each vertical side of a cell may contain at most two diagonals. 
The horizontal visibility decomposition of F, denoted by HD{F), is defined similarly by extending 
the horizontal edges of V. 

Our v-v-map is on VD[F), i.e., for each cell C of VD[F), the link distances of the min-link v-v- 
paths from s to all points in C that are not on diagonals are the same and we denote this distance 
by diSyy{C), and for each diagonal d of VD{F), the link distances of the min-link v-v-paths from 
s to all points on d are the same and we denote this distance by disyy{d). In fact, it holds that 
disyy{d) = min{diSyy{Ci),disyy{Cr)}, where Ci and Cy are the two cells on the left and right of d, 
respectively. The goal of our algorithm for computing the v-v-map is to compute diSyy{C) for each 
cell C and dis{d) for each diagonal d of VD{F). In addition, we also need to maintain some path 
information to retrieve an actual path for each query. 

Similarly, our h-v-map is also on VD(F), but the h-h-map and the v-h-map are both on HD[F). 

In Section [2Tl we review and discuss the DN algorithm in a way that will help us to introduce 
our algorithm later. The DN algorithm only labels the diagonals of VD[F). Das and Narasimhan [9] 
claimed that queries can be answered by using only the distance values on the diagonals. However, 
this is not clear to us. We, instead, also label the cells as discussed above. We will show in Section 
12.21 that the DN algorithm can be easily adapted to computing the link distances for the cells too. 
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Fig. 2. Illustrating the vertical visibility decomposition of T . 

We will also discuss in Section 12.21 how to maintain path information to retrieve an actual path for 
each query. 

2.1 The DN Algorithm 

The DN algorithm also computes the four maps discussed above. We consider the v-v-map first. 
The algorithm works on the vertical visibility decomposition VD{T') and the goal is to compute 
disyy{d) for each diagonal d. To simplify the notation, unless otherwise stated, we use dis{-) to refer 
to disyy{-). Initially, all diagonals have distance value oo except dis{ds) = 1, where dg is the diagonal 
through s. Note that if a diagonal d is on a side e of a cell, then whenever dis{d) is updated, dis{e) 
is automatically set to dis{d). 

The DN algorithm has many phases. In the i-th phase for i > 0, the algorithm determines 
the set Vi of diagonals d whose distances dis{d) are equal to 2 i + 1, and these diagonals are then 
“labeled” with distance 2z +1 (e.g., see Fig. [3]). Initially, i = 0, and Vq consists of the only diagonal 
through s. As discussed in [9], if we put light sources on the diagonals in Vi_i, then Vi consists of 
all new diagonals that will get illuminated with light emanating horizontally from the light sources. 

Consider a general z-th phase for z > 1. We assume Vi-i has been determined. There are two 
procedures: right-sweep and left-sweep. In the right-sweep (resp., left-sweep), we illuminate the 
diagonals in the direction from left to right (resp., from right to left). The right-sweep procedure 
starts from the locally-rightmost diagonals of I^_i, dehned as follows. Consider any diagonal d in 
Vi-i. Let C be the cell of VD{F) on the right of d, i.e., d is on the left side of C. Let be the 
right side of C. If dis{er) 7 ^ 2z — lo, then d is a locally-rightmost diagonal of Vi-i. Similarly, the left- 
sweep starts from the locally-leftmost diagonals of Vi-i. Both locally-leftmost and locally-rightmost 
diagonals are referred to as locally-outmost diagonals. Below, we first discuss the right-sweep. 

For each locally-rightmost diagonal d, we put a rightward “light beam” on d, and let B{d) 
denote the beam; we also call d the generator of the beam. Formally, we may dehne the beam 
as the segment d associated with the rightward direction. Initially we insert all locally-rightmost 
diagonals into a min-heap whose “keys” are the x-coordinates of the diagonals (i.e., the leftmost 
diagonal is at the root). By using Hr, the diagonals involved in the right-sweep will be processed 
from left to right. As discussed in [9], this can control the right-sweep in a global manner and thus 
avoid the same diagonal being processed many times. Although initially each diagonal of Hr has 

^ Alternatively, one may replace dis{er) 7 ^ 2 i — 1 by dis{er) > 2 i — 1 for the definition, but we choose to use 
dis{er) 7 ^ 2 z — 1 only for making our discussion later easier and this will not affect the algorithm running time. 
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Fig. 3. Illustrating the DN algorithm: all diagonals have been labeled. Fig. 4. Illustrating a cell C where d is on the 
Note that the three thick dash-dotted diagonals (labeled 5) are swept left side of C. 
twice in the second phase. 

only one beam, later we will insert more diagonals into Hu and a diagonal d may have more than 
one beam, in which case B[d) denotes the set of beams on d. 

As long as Hr is not empty, we repeatedly do the following. 

We obtain the leftmost diagonal d of Hr and remove it from Hr. Let C be the cell on the right 
of d (e.g., see Fig. [3]). We process d in the following way. Intuitively we want to propagate the beams 
of B{d) to other diagonals in C. Let ei and denote the left and right sides of C, respectively. 
Note that d is on e;. Recall that each cell side has at most two diagonals. 

If ei has another diagonal d (e.g., see Fig. [3|) and d has not been labeled (i.e., dis{d) = oo), 
then we set dis{d) = 2i + 1. The beam set B{d) of d is set to 0 since there is no beam from B(d) 
illuminating d. Further, although B(d) = 0, we associate the leftward direction with it, because 
if later the left-sweep procedure does not illuminate d, d will be a locally-leftmost diagonal of Vi 
and generate a leftward beam in the next phase. We “temporarily” mark d as a locally-leftmost 
diagonal. We say “temporarily” because d may not be locally-leftmost any more after the left-sweep 
procedure, as discussed later. 

If d has been labeled, then dis{d) must be 2i + 1. To see this, because we are currently at the 
Lth phase, any diagonal that has been labeled must have a distance value at most 2i + 1. But if 
dis{d) <2i + l, then d has already been processed, and since d and d are on the same cell side, we 
must have already obtained dis{d) = dis{d) < 2i + 1, contradicting with dis{d) = 2i + 1. Hence, 
dis{d) = 2i + 1. In this case, we do nothing on d. 

Next, we consider the diagonals on the right side of C. Depending on the values of dis{er), 
there are several cases. Since we are at the i-th phase, either dis{er) = oo or dis{er) <2i + l. 


1. If dis{er) < 2i + 1, then we do nothing on Cr because all diagonals on have already been 
processed. 

2. If dis{er) = oo, then we set dis{er) = 2z -|- 1. If Cr does not have any diagonals, we are 
done. Otherwise, for each diagonal d' on e^, we determine the portions of beams of B{d) that 
can illuminate d', which are the rightward projections of B{d) on d' (beams of B{d) may be 
“narrowed” or “split”, e.g., see Fig. [5]). We use B{d) D d' to denote the above portions of B{d). 
It is possible that B (d) fl d' is empty. 


6 























































(a) (b) (c) 


Fig. 5. Illustrating some beam operations in a right-sweep procedure: In (a) and (b), beams are split, and some beams 
are “narrowed” and some beams “terminate” at obstacle edges; in (c), beams are merged. 

If B{d) n d' = 0, then we temporarily mark d! as a locally-rightmost diagonal and set B{d') = 0 
with the rightward direction] otherwise, we set B{d') = B{d) d' and insert d' to Hu (later we 
will propagate the beams of B{d') further to the right). 

3. If dis{er) = 2i-|-l, this case happens because was illuminated by beams from another diagonal 
d on e^. Hence, each diagonal d' on Cr may already have a non-empty B[d'). But d' may receive 
more beams from B (d). We first determine B (d) fl d' (as defined above) and then do a “merge” 
operation by merging B(d) D d' with B{d') (e.g., see Fig. [5]). Finally, we set B{d') to the above 
merged set of beams (with the rightward direction). Due to the merge operation, B{d') may 
contain more than one beam. 

If B{d') was empty before the merge and becomes non-empty after the merge, then we insert 
d' into Hpi. If B{d') is still empty after the merge, then we temporarily mark d' as a locally- 
rightmost diagonal. If B[d') was non-empty before the merge, then d' is already in Hu, so we 
do not need to insert it into Hu again. 

The above finishes the processing of d. The right-sweep is done once Hr becomes empty. For 
the implementation, we need to maintain the beams of B(d) on each diagonal d. To this end, as in 
[9], we can use a balanced binary search tree (e.g., a 2-3-4 tree [8]) such that we can perform the 
“merge”, “split”, “insertion”, “deletion”, and “search” operations each in logarithmic time. 

The left-sweep procedure is similar. However, there is one subtle thing. If the sweep encounters 
a diagonal d that has been labeled by the right-sweep, then this is ignored and we keep sweeping 
as if d were not labeled. As discussed in [9], the reason for this is that the left-sweep may reach 
more cells than the right-sweep (e.g., see Fig. [3]). In this way, each diagonal can be processed at 
most twice in a phase. But no diagonal can be processed in more than one phase. Also, suppose 
a diagonal d was temporarily marked as a locally-outmost diagonal during the right-sweep; if d 
is illuminated again in the left-sweep but d is not marked locally-outmost in the left-sweep, then 
we clear the previous mark on d (i.e., d is not considered locally-outmost any more). After the 
left-sweep, the remaining locally-outmost diagonals are exactly those that will be used in the next 
phase. 

The above describes the f-th phase of the algorithm. The algorithm is done after all diagonals 
are labeled. Due to that heaps are used to control the sweep procedures and balanced binary search 
tree are used to support beam operations, the total time of the algorithm is 0(n log n) because each 
diagonal can be processed at most twice in the entire algorithm (once by a left-sweep procedure 
and once by a right-sweep procedure). Clearly, the space is 0{n). 

For computing the h-v-map, the algorithm is similar except that Vq now consists of all diagonals 
that intersect d(,, where d'g is the horizontal line segment extended from s until it hits the obstacles. 
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and in each i-th phase for i > 0, the value dishvid) is set to 2{i + 1) for any diagonal d in Vi. The 
algorithms for computing h-h-map and v-h-map in HD{T) are symmetric. 

2.2 Labeling Cells and Retrieving Paths 

We first show how to slightly modify the DN algorithm to label cells. We only discuss the v-v-map 
since the other maps are similar. 

Recall that for each cell C of all the points in C that are not on the diagonals have the 

same distance values. Hence, to compute diSwiC), it is sufficient to know the distance value for 
any arbitrary point in C that is not on any diagonal. To this end, for each cell C, we add a vertical 
segment in C with its upper endpoint on the upper side of C and its lower endpoint on the lower 
side of C such that the segment is not overlapped with the left side or the right side of C; hence, 
no point of the segment is on any diagonal, and we call the segment a “fake diagonal”. Let VD'{T) 
denote VD[T) with all fake diagonals. We run the DN algorithm on VD'[F) and treat all fake 
diagonals as “ordinary diagonals” to label all diagonals and fake diagonals of VD'{F). Finally, we 
label each cell of VD{T) with the same distance value as its fake diagonal. Since the total number 
of fake diagonals are 0{n), the running time is still 0(nlogn) and the space is 0(n). 

To obtain an actual min-link v-v-path from s to t, we need to maintain additional information 
on our v-v-map. No details on this are given in [9]. We briefly discuss it below for the completeness 
of this paper. Essentially, when we label the cell sides (and the fake diagonals), we need to record 
how we reach there. Specifically, suppose we label a diagonal d on the right side of a cell in a 
right-sweep procedure due to a beam from a diagonal on the left side of the cell; then we record 
any such beam at d (it is sufficient to record any point on d in the beam) along with its generator, 
and in the case where B{d) is empty, d is a locally-outmost diagonal and the path should make a 
turn there. Further, for each locally-outmost diagonal d, it is reached by a beam that illuminates 
the cell side containing d and we record that beam for d so that we know the turn on d is for that 
beam. With this path information, for any query point t, we can easily retrieve an actual min-link 
v-v-path from t back to s in time we claimed before. We omit the details. 

3 The Main Scheme of Our Algorithm 

In this section, we focus on the main scheme of our algorithm while leaving the algorithms for 
processing corridors in Section 01 We first introduce the corridor structure. 

3.1 The Corridor Structure 

The corridor structure in rectilinear domains is similar to that in general polygonal domains m- 
Let be the dual graph of the vertical visibility decomposition VD{F) (see Fig.[6]), i.e., each node 
of Gytd corresponds to a cell of VD{F), and each edge of Gytd connects two nodes corresponding 
to two cells sharing the same diagonal. Based on Gytd^ we obtain a corridor graph G as follows 
(see Fig. ED. First, we remove every degree-one node from Gytd along with its incident edge; repeat 
this process until no degree-one node exists. Second, remove every degree-two node from Gytd and 
replace its two incident edges by a single edge; repeat this process until no degree-two node exists. 
The remaining graph is G. The cells in VD{T) corresponding to the nodes in G are called junction 
cells (see Fig. ED- We consider the diagonal through s as a degenerate junction cell. Similar to 
the corridor structure in the general polygonal domains [20], the graph G has 0{h) nodes and 




Fig. 6 . Illustrating the vertical visibility decomposition 
(the dashed segments are diagonals) and its dual graph 
Gvtd- 


Fig. 7. Illustrating the graph G, and the corridor (shaded 
by slashes) bounded by Pi and Pj. 


0(h) edges. The removal of all junction cells from VD(T) results in 0(h) corridors, each of which 
corresponds to an edge of G. 

The boundary of any corridor C consists of four parts (see Fig. [7]): (1) The boundary portion 
of an obstacle Pi, from a point a to a point 6; (2) a diagonal be; (3) the boundary portion of an 
obstacle Pj from c to a point d; (4) a diagonal da. The two diagonals be and ad are called the doors 
of the corridor C. The corridor C is a simple rectilinear polygon. 

3.2 The Main Idea 

We focus on computing the v-v-map on VD(P), and the other three maps can be computed similarly. 
Our goal is to label d for each diagonal d of VD(P), i.e., compute the distance value dis(d) = 
disyy(d). As discussed in Section [2j2l the same algorithm can be used to label cells as well. As in 
the DN algorithm each diagonal d will maintain a beam set B(d). 

Here is some intuition that motivates us to improve the DN algorithm. When we are running 
the DN algorithm, a sweep procedure will enter each corridor though one of its two doors, and 
the procedure will either sweep the entire corridor and leave the corridor through the other door, 
or terminate inside the corridor (in which case the sweep “hits” another sweep that entered the 
corridor through the other door and both sweeps terminate after “collision”). This means that if 
we can determine the beams and the distance values at the doors of a corridor, then we can process 
the corridor independently in a more efficient way since the corridor is a simple rectilinear polygon. 

Our algorithm follows the similar scheme as the DN algorithm. The sweep in the junction cells is 
still processed and controlled in a global manner in each phase. However, whenever the sweep enters 
a corridor through one of its doors, the corridor will be processed independently by using our more 
efficient corridor-processing algorithm given in Section [4] (one may view that the sweep procedure 
jumps from one junction cell to another through the corridor connecting them). We should point 
out that since in the DN algorithm a diagonal may be processed twice in the two sweep procedures 
in the same phase, an entire corridor may be processed twice in the same phase (that case happens 
only if the beams on a door can illuminate the other door directly, and vice versa). 

The running time of our algorithm is 0(n + hlogh). More specifically, since there are 0(h) 
junction cells, the time spent on processing the diagonals in all junction cells is 0(hlogh), and the 
processing on all corridors takes 0(n) time because the number of vertices of all corridors is 0(n), 
in addition to another 0(/i log h) time spent on maintaining the beams on all diagonals. 
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3.3 The Algorithm 


To help the reader understand the algorithm, we first describe the hrst few steps of the algorithm 
and then delve into the details of the general steps. 

Initially, we set dis{d) to oo and B{d) = 9 for each diagonal d except that dis{ds) = 1 and 
B(ds) = {ds}) where dg is the diagonal through s. In our discussion of the corridor structure, dg is 
considered as a degenerate junction cell. Here, for ease of discussion, we assume dg is in the interior 
of a corridor C(s) as shown in Fig. [71 We first process C(s), i.e., label all diagonals in C(s). For 
the purpose of describing our algorithm, one may assume we still use the DN algorithm to process 
C(s), and later we will replace the DN algorithm by our new and more efficient corridor-processing 
algorithm in Section 01 

Denote by and Dy the two doors of C(s) (e.g., see Fig. [7]). Suppose after the above processing 
of C(s), dis{Dx) = 2 i -|- 1 and dis{Dy) = 2j + 1 for some integers i and j. Since dis{Dx) is obtained 
“locally” in C{x), i.e., dis{Dx) is the link distance of the local min-link v-v-path from s to Dx 
in C(s), it may not be “set correctly”, i.e., it may not be the link distance of the global min-link 
v-v-path from s to Dx in B. The value dis{Dy) has the same issue. However, we show below the 
at least one of dis{Dx) and dis{Dy) must have been set correctly. 

Without loss of generality, assume j < i. 

Observation 1 If j <i, then the value dis{Dy) has been set correctly. 

Proof. Let vr be any global min-link v-v-path from s to Dy in D. For any subpath tt' of vr, we use 
dis{'K') to denote the link distance of tt'. If vr is in C(s), then our observation follows. Otherwise, 
vr must cross one of the doors of C(s). Let p be the first point on a door of the corridor if we go 
from s to Dy along tt, and let 7r{s,p) denote the sub-path of vr from s to p. If p is on Dy, then 
7r(s,p) is a min-link v-v-path from s to Dy. Since 7r(s,p) G C(s) and dis{Dd) is the link distance 
of the local min-link v-v-path from s to Dy in C(s), we have disfn) > dis{'K{s,p)) > dis{Dy), 
and thus our observation follows. If p is on Dx, then dis{7r(s,p)) = 2i + 1 > 2j + 1. Therefore, 
dis{'K) > dis{'K{s,p)) >2j + l = dis{Dy) and our observation also follows. □ 

Another way to see the observation is that if we were running the DN algorithm, after the 
(j — l)-th phase, neither Dx nor Dy is labeled, and no diagonal outside the corridor is labeled 
either. In the j-th phase, Dy will be labeled and thus dis{Dy) = 2j + 1 is set correctly. 

After the processing of C(s), the beams on Dy, i.e., B(Dy), have also been obtained. The next 
step is to process the diagonal Dy by propagating the beams of Dy outside the corridor. 

Next, we describe the details of the general steps of our algorithm. 

We use a min-heap H to maintain the diagonals in all junction cells where the “keys” are the 
distance values the diagonals currently have (and these values may not be set correctly), with the 
smallest key at the root of H. Since there are 0{h) junction cells, the size of H is always 0{h). 
Each diagonal d in H is also associated with its beam set B{d) (along with its direction). As in 
the DN algorithm, it is possible that B{d) is empty, in which case d might be a locally-outmost 
diagonal, but it is also associated with a direction for generating a beam towards that direction in 
the next phase. 

If some diagonals of H have the same keys, we break the ties by applying the following rules. 
Consider two diagonals di and d 2 in H with dis{di) = dis{d 2 ). If B{di) and B{d 2 ) are both empty 
or both non-empty, then we break ties arbitrarily. Otherwise, without loss of generality, assume 
B{di) is not empty but B{d 2 ) is empty. Then, we consider the key of di to be smaller than that 
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of d 2 in H. The reason is as follows. Since B{di) and B{d 2 ) = 0, the current sweep procedure 
should be over before processing d 2 while the sweep should continue after processing di, and thus, 
we should process di before d 2 - Therefore, our way of resolving ties in H is crucial and consistent 
with the DN algorithm. In the following, to differentiate from keys of other heaps, for any diagonal 
d, even if d is not in H, we consider dis{d) along with B{d) as the global-key of d, and whenever 
we compare the global-keys of diagonals, we always follow the above rules. 

Consider any corridor C with two doors d and d'. Suppose the beams of B{d) are going inside 
C, and we want to process C (i.e., compute the v-v-map in C) using the beams of B{d). We say the 
above way of processing C is in the direction from d to d'. As will be seen later, a corridor may be 
processed twice: once from d to d' and the other from d' to d. Due to the special geometric structure 
of the corridor, we have the following observation that will be useful later. 

Observation 2 Suppose d and d' are the two doors of a corridor C, and the direction of the beams 
of B{d) is towards the inside ofC. Then after C is processed by using B{d), the beam set of d' is 
not empty. 

Proof. Let VD{C) denote the vertical visibility decomposition of the corridor C. Consider the cell 
C of VD{C) that contains d'. Without loss of generality, assume d' is on the right side of C. Denote 
by Cr the right side of C. 

We claim that d' is the entire right side of C, i.e., d' is e,.. Indeed, according to our corridor 
structure, d' is an extension of a vertical obstacle edge e and one endpoint p of e is also an endpoint 
of d' and the other endpoint of e is outside the corridor C. This means p is also an endpoint of 
Cr and the other endpoint q of d' than p is not on e but on another obstacle edge e'. Due to our 
general position assumption that no two vertical edges are collinear, q must be in the interior of e', 
implying that q is also an endpoint of e^. Hence, we obtain d' = Cr = pq. 

Note that we obtain the beam set of d' from the rightward beams of the diagonals on the left 
side of C. Now that d' is the entire right side of C, d' will receive all beams of any diagonal on the 
left side of C. Hence, the beam set of d' cannot be empty. □ 

To show the correctness of our algorithm, we will argue that our algorithm is consistent with 
the DN algorithm. We will show that after the algorithm finishes, each diagonal in any junction 
cell is correctly labeled, i.e., both its distance value and its beam set are the same as those in the 
DN algorithm. Let d* be the diagonal in the root of H. Our algorithm will maintain the following 
three invariants. 

1. The diagonal d* is correctly labeled. Further, for any other diagonal d in a junction cell, if the 
global-key of d is no larger than that of d*, then d has been correctly labeled. 

2. For any diagonal d in a junction cell, if dis{d) oo and the global-key of d is larger than that 
of d*, then d is in H. 

3. For any corridor C with two doors d and d', if C is processed in the direction from d to d', then C 
will never be processed from d to d' again in the entire algorithm (although C may be processed 
later in the other direction from d' to d). 

Initially H = ^. Recall that dis{ds) = 1 and B{d) = {d^}. We consider the diagonal dg through 
s as a degenerate junction cell. Specifically, we consider dg as two duplicate diagonals with one 
generating a rightward beam and the other generating a leftward beam from the entire dg. We 
insert these two diagonals into H. Clearly all algorithm invariants hold. In the following we will 
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describe the details of our algorithm. To avoid the tedious discussion, we will not explicitly explain 
that the algorithm invariants are maintained after each step, but rather discuss it briefly later after 
the algorithm description. 

As long as H is not empty, we repeatedly do the following. 

By an extract-min operation on H, we obtain the diagonal d* of H with the smallest global-key 
and remove it from H. We assume dis{d*) = 2i -\-l for some integer i. If we were running the DN 
algorithm, we are currently working on the i-th phase. There are two cases depending on whether 
B{d*) = 0. 

3.3.1 B{d*) ^ 0 

We first discuss our algorithm for the case where B{d*) ^ 0. If we were running the DN algorithm, 
we would be in the “middle” of the i-th phase (because B{d*) implies that all locally-outmost 
diagonals have already been processed). So we should continue the left-sweep and right-sweep of the 
i-th phase. The first question is where we should start the sweep. Let S be the set of all diagonals 
in H that have the same global-keys as d* (according to our way of comparing global-keys, for each 
diagonal d G S', it holds that dis{d) = dis{d*) and B{d) ^ 0). By the first algorithm invariant, 
all diagonals in S have been correctly labeled. Due to our corridor structure, the sweeps of the 
i-th phase “paused” at the diagonals in S. To continue the i-th phase, we “resume” the sweeps 
from these diagonals. Unlike the DN algorithm where we complete the right-sweep before we start 
the left-sweep, here, before the pause, we may have already done some left-sweep and right-sweep. 
Hence, the two sweeps may be somehow “interleaved” and our algorithm will need to take care of 
this situation. 

Due to our way of breaking ties in H, the diagonals of S can be found by continuing the 
extract-min operations on H in 0(151 log |LI|) time (i.e., all diagonals of S are removed from H). 
We also let S contain d*. Let Sr (resp., Sr) be the subset of the diagonals of S whose beams are 
rightward (resp., leftward). Intuitively, the right-sweep (resp., left-sweep) paused at the diagonals 
in Sr (resp., Sr), and thus, we resume it from the diagonals in Sr (resp., Sr). Below we focus on 
the right-sweep, and the left-sweep is very similar. 

Recall that in the right-sweep of the DN algorithm we use a heap Hr to guide the procedure. 
Here we do the same thing and process the diagonals of Hr from left to right. We build a heap Hr 
by inserting the diagonals of Sr, and the “keys” of diagonals in Hr are their x-coordinates such 
that the leftmost diagonal is at the root. (Similarly, we build a heap Hr on Sr for the left-sweep.) 

The algorithm essentially performs the i-th phase as the DN algorithm. But since here the right- 
sweep and left-sweep may be interleaved, in the following discussion, some diagonals may have two 
sets of beams with opposite directions. However, our algorithm makes sure that if a diagonal d has 
two sets of beams with opposite directions, it will not be in H (i.e., it has been removed from H), 
but in both Hr and Hr d has not been processed yet. To differentiate the two sets of beams, we 
use Br{d) (resp., Bi{d)) to denote the beam set of any diagonal d in Hr (resp.. Hr), meaning that 
the direction of the beams is rightward (resp., leftward). 

During the right-sweep, if we find a new diagonal d that has the same global-key as d* , then d 
will be inserted to Hr and d will be removed from H if it is already in H. Hence, all diagonals of 
Hr have the same global-key as d* . In fact. Hr maintains all diagonals in junction cells that will 
be processed in the subsequent right-sweep of the i-th phase. 

As long as Hr is not empty, we repeatedly do the following. 
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We obtain the leftmost diagonal d of Hji (which is at the root of Hji) and remove it from Hji. 
The beams of Br{d) may enter either a junction cell or a corridor. In the sequel, we discuss how to 
process d in these two cases. 

The beams of Br{d) entering a junction cell. Let C denote the junction cell that the beams of Br{d) 
enter. In this case, our way of processing d is similar to the DN algorithm. Let ei and denote the 
left and right sides of C, respectively. Note that d is on the left side e/. Recall that each cell side 
may have two diagonals. To process d, we update the labels of all other diagonals of C, as follows. 
Suppose there is another diagonal d on ei (e.g., see Fig. SD- 

If dis{d) = oo, we set dis{d) = dis{d) and B{d) = 0 with the leftward direction. Then, we insert 
d into H. 

If dis{d) 7 ^ oo but dis{d) > dis{d), then we set dis{d) = dis{d) and B{d) = 0 (with the leftward 
direction). Since dis{d) was greater than dis{d) but not oo, by the second algorithm invariant, d 
is already in the heap H. Hence, after dis{d) and B{d) are reset as above, we do a “decrease-key” 
operation on d in LI since the global-key of d has been decreased. 

If dis{d) < dis{d), we do nothing. 

Next, we consider the diagonals on the right side Cr of C. Depending on the distance value 
dis{er), there are several cases. Note that if dis{er) 7 ^ 00 , then automatically got the value 
dis{er) because a diagonal d' on was labeled with the same distance value. 

1 . If dis{er) = 00 , then we set dis{er) = 2i -|-1. If does not have any diagonals, then we are done 
with processing d. Otherwise, for each diagonal d' on e^, we set dis{d') = 2i + 1 and determine 
Br{d) n d', i.e., the portion of Br{d) that can illuminate d'. 

If Br{d) n d' 7 ^ 0, then we set Br{d') = Br{d) fl d' and insert d' to Hr (note that d' has the same 
global-key as d*). Otherwise, we set B{d') = 0 with the rightward direction and insert d' into 
H, because the global-key of d' is larger than that of d* and the right-sweep should stop at d'. 
Note that before the insertion, d' was not in H as dis{d') was 00 . 

2. If dis{er) 7 ^ 00 but dis{er) > 2i + 1, then the algorithm is similar as above. For each diagonal 
d' on Cr, dis{d') was equal to dis{er), and now we set dis{d') = 2i-\-l and determine Br{d) D d'. 
If Br{d) n d' 7 ! 0, we set Br{d') = Br{d) fl d' and insert d' into Hr. Since dis{d') was equal to 
dis{er), the global-key of d' was larger than that of d* and dis{d') was not 00 . By the second 
algorithm invariant, d' was already in H. Hence, we remove d' from H. Note that we can perform 
a remove operation on d' in H by doing a decrease-key operation followed by an extract-min 
operation [ 8 ]. 

If Br{d) n d' = 0, we set B{d') = 0. Since d' was already in H, after setting dis{d') = 2i + l and 
B{d') = 0, we need to do a decrease-key operation on d' in H. 

3. If dis{er) < 2i -|- 1, we do nothing. 

4. If dis{er) = 2i + l, due to that the left and right sweeps are interleaved, may have got labeled 
from the right-sweep, the left-sweep, or both. We discuss the three cases below. The algorithm 
for this case is also simple, but we need a few more words to explain why the algorithm works 
in that way. 

(a) If Cr got the value dis{er) by the right-sweep only, as in our discussions in the right-sweep 
of the DN algorithm, this case happens because was illuminated by beams from another 
diagonal d on the left side ei. Hence, for each diagonal d' on e^, dis{d') and Br{d') have been 
set. Further, if Bj.{d') = 0, then d' is in id; otherwise d' is in Hr. 
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For each diagonal d' on e^, we first determine Br{d)rid' , and then set Br{d') by merging the 
original Br{d') with Br{d) n d' . If Br{d') was non-empty before the merge, we do nothing 
since it is already in Hji. If Bj.{d') was empty before the merge and becomes non-empty 
after the merge, then we insert d' into Hji and remove it from H. If Br{d') is still empty 
after the merge, then we do nothing since it is already in H. 

(b) If Br got the value dis{er) by the left-sweep only, then was illuminated by beams from its 
right. As discussed in the DN algorithm, if the right-sweep procedure sweeps a diagonal d' 
on Br, the fact that d' was swept already by the left-sweep procedure should be ignored in 
the sense that we should keep propagating the beams of Br{d) to the right of d'. The details 
are given below. 

For each diagonal d' on e^, we first determine Br{d) Cid'. 

— If Br{d) n d' 7 ^ 0, then we set Br{d') = Br{d) fl d' and insert d' into Hr. If d' was already 
in H, then we remove it from H. 

— If Br{d) n d' = (Ji, no beam from Br{d) illuminates d'. At first sight, it seems that we 
should insert d' into H with an empty beam set and the rightward direction. Below we 
elaborate on whether we should do so. 

Note that d' is a door of a corridor C that is locally on the right of d'. It is possible that 
Br got labeled because the left-sweep was from C, in which case Bi{d') is not empty by 
Observation [21 and thus, d' must be already in Hi (since dis{d') = 2i + 1) and we do 
not need to insert d' to H because d' will be processed in the left-sweep of the current 
phase. 

On the other hand, if b^ has another diagonal d” , then it is possible that got labeled 
because of the processing of the corridor on the right of d" during the left-sweep. In this 
case, d' may have got labeled because of the processing of d", in which case as in the DN 
algorithm the beam set of d' must be empty and rightward, and thus d' is already in H 
with B{d') = 0 and the rightward direction and we do not insert d' into H again. But if 
d' has not been labeled yet, then we insert d' into H with B[d') = 0 and the rightward 
direction. Therefore, for the case where has another diagonal d” , if d' is already in H, 
we do nothing; otherwise we insert d' to H. 

It is also possible that Br got labeled “simultaneously” because of the processing of the 
two corridors on the right of d' and d” , in which case by Observation [2] we again have 
Bi{d') 7 ^ 0, and thus d' is already in Hi. Hence, we do not need to insert d' to H. 

In summary, for the case Br{d) D = 0, if d' is in neither Hi nor H, then we insert d' 
into H with B[d') = 0 and the rightward direction; otherwise we do nothing. 

(c) If Br got the value dis{Br) by both the right-sweep and the left-sweep, this is a combination 
case of the above two cases. 

For each diagonal d' on e,., we determine Br{d) D d' , and then set Br{d') by merging the 
original Br{d') with Br{d) fid'. 

If Br{d') was non-empty before the merge, then we do nothing since d' is already in Hr. 

If Br{d') was empty before the merge and becomes non-empty after the merge, then we 
insert d' into Hr. Further, if d' is in H, then we remove it from H. 

If Br{d') is still empty after the merge, as in the above second case, we do the following. 
If d' is in neither Hi nor dd, then we insert d' into H with B{d') = 0 and the rightward 
direction; otherwise we do nothing. 

We are done with processing d when the beams Br{d) of d enter a junction cell. 
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The beams of Br{d) entering a corridor. Let C denote the corridor that the beams enter. We process 
C using the beams of Br{d). Again, one may assume we still use the DN algorithm to process C, 
and later we will replace it by our corridor-processing algorithm in Section 01 

Let 6 be the distance value labeled on the other door d' of C by the above processing and let B' 
denote the corresponding beam set on d'. Let dis{d') and B{d') be the original distance value and 
beam set at d' before the processing of C. By the third algorithm invariant, this is the first time 
C is processed in the direction from d to d'. Hence, if dis{d') ^ oo, then the value dis{d') must be 
obtained by the sweep from outside C, i.e., beams in B{d') are towards the inside of C. 

Due to the above processing of C, we have obtained another distance value 6 and beam set B' 
for d' . Hence, we need to update the label of d' and possibly insert d' to some heap. Depending on 
the value of dis{d'), there are several cases. 

1. If dis{d') is oo, then we set dis{d') = 5 and B{d') = B' . If 5 > 2i -|- 1, then we insert d' into 

H.li 5 = 2i + 1, since dis{d) = 2i + 1, d' must be illuminated directly by the beams in Br{d) 
and the beams of B[d') are still towards right. By Observation [2l B' ^ 0. Hence we obtain 
Br{d') = 7 ^ 0 (we set Br{d') to B' because the beams of B' are rightward). Finally, we insert 

d' into Hr. 

2. If dis{d') <2i + \, then the global-key of d' is smaller than that of d* because dis{d*) = 2i + l. 
By the first algorithm invariant, d' has been correctly labeled. Recall that the direction of B[d') 
is towards the inside of C. Also by the first algorithm invariant, d is correctly labeled, and the 
direction of the beam set of Br{d) is towards the inside of C. This means we have computed 
complete information on the two doors of C for the min-link v-v-paths from s to the points 
inside C. Hence, we can do a “post-processing” step to compute the v-v-map in C by using the 
beams of Br{d) and B[d'). We will give a corridor-post-processing algorithm for this step later 
in Section m 

3. If dis{d') = 2i -\- 1, then d' has been labeled in the current phase. 

If B[d') 7 ^ 0, then the global-key of d' is the same as that of d* . By the first algorithm invariant, 
d' has been correctly labeled. Then, as above, since both d and d' have been correctly labeled, 
we do a “post-processing” to compute the v-v-map in C using Br{d) and B{d'). 

If B{d') = 0, then the global-key of d' is strictly larger than that of d*. By the second algorithm 
invariant, d' is already in H. 5 > 2i -\- 1 , then we do nothing. If <5 = 2i -|- 1, then as in the 
above first case, we set Bj.{d') = B' ^ finally, we insert d' into Hr and remove d' from H. 

4. The remaining case is when dis{d') 7 ^ 00 and dis{d') > 2z-|- 1. We claim that this case can never 
happen. 

Indeed, assume to the contrary that this case happens. Recall that the beams of B{d') are 
towards the inside of the corridor C. Let C be the junction cell that contains d'. Without loss 
of generality, assume d' is on the left side of C. Clearly, d' got labeled after some diagonal d” in 
C was processed. Since the beams of d" that illuminate d' must be towards the cell C, they are 
from the corridor that is bounded by d". By Observation [21 the beam set of d" is not empty. 
Hence, it must be the case that dis{d") = dis{d') > 2i + 1. Since we use the heap H to guide 
the main algorithm and we are currently processing the diagonal d with dis(d) = 2i -\- 1, all 
diagonals of junction cells that have been processed must have distance values at most 2i + 1 . 
However, the above shows that the diagonal d" has been processed with dis{d") > 2i -|- 1, 
incurring contradiction. Therefore, the case where dis{d') 7 ^ 00 and dis{d') > 2i -|- 1 cannot 
happen. 
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The above finishes the processing of the diagonal d. The right-sweep procedure is done after 
the heap Hfi becomes empty. Afterwards we do the left-sweep from the diagonals of 5^ using the 
heap Hl in the symmetric way. We omit the details. 

This finishes our algorithm in the case in which B(d*) is not empty for the diagonal d* at the 
root of H. We briefly discuss that after the above processing our algorithm invariants still hold. 

Indeed, before the diagonals of Sji are processed, if we were running the DN algorithm, at some 
moment during the f-th phase, the diagonals of Sr would be labeled the same as in our algorithm. 
Our algorithm processes the diagonals of Sr in a way consistent with the DN algorithm. Based on 
our previous discussion on how the sweep procedures of the DN algorithm can sweep the corridors, 
after all diagonals of Sr are processed, the diagonals in H with the smallest global-key must have 
been correctly labeled because if we were running the DN algorithm, at some moment the algorithm 
would label them in the same way. Also, since we process diagonals in the order of their global-keys, 
any diagonal that has global-key smaller than that of the root of H must have been processed and 
labeled correctly. Hence, the first algorithm invariant follows. For the second invariant, whenever 
a diagonal has its distance value set to non-infinity for the first time, it is always inserted into H. 
For the third invariant, as discussed in the algorithm description, a corridor C is processed only if 
a door d of C is processed and beams of B(d) are towards C. Although d may be processed in the 
algorithm twice, it is processed only once in either the right-sweep or the left-sweep. Hence, with 
beams towards the inside of C, d is only processed once in the entire algorithm, implying that C is 
processed only once in the entire algorithm in the direction from d to the other door of C. 


3.3.2 B{d*) = 0 

In the sequel, we discuss the case where B(d*) = 0. The algorithm is simpler in this case. First, 
we find the set S of diagonals in H that have the same global-key as d*, which can be done 
in 0(|5| log |id|) time by keeping doing the extract-min operations (i.e., all diagonals of S are 
removed from H). We also let S contain d*. According to our way of comparing global-keys, all 
diagonals of S have empty beam sets. If we were running the DN algorithm, the diagonals of S 
would be locally-outmost and we would be about to start the {i + l)-th phase (not the i-th phase). 
As in the previous case, we run the two sweep procedures starting from the diagonals of S. Let Sr 
(resp., Sr) be the subset of diagonals of S whose beam directions are rightward (resp., leftward). 
We build a min-heap Hr (resp.. Hr) on the diagonals of Sr (resp., Sr). Below, we only discuss the 
right-sweep since the left-sweep is similar. 

Since we are doing the right-sweep in the {i + l)-th phase, each diagonal of Sr will generate 
a beam from the entire diagonal, and all new diagonals illuminated in the right-sweep will get 
distance value 2i + 3 instead of 2i -|- 1. From now on, we associate each diagonal of Sr with the 
beam, i.e., for each d € Sr, Br{d) now consists of the only beam on d although it was empty in 
the Tth phase. Since each diagonal of Sr originally got an empty beam set (at the end of the Tth 
phase), by Observation O the beam of Br{d) cannot be towards a junction cell and thus it must be 
towards the inside of a corridor. 

As long as Hr is not empty, we repeatedly do the following. 

We obtain the leftmost diagonal d of Hr (which is at the root) and remove it from Hr. Let C 
denote the corridor that the beams of Br{d) enter. We process the corridor C using the beams of 
Bj.{d). Again, one may assume we still use the DN algorithm to process C, and later we will replace 
it by our corridor-processing algorithm in Section HI 
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Let 6 be the distance on the other door d' obtained by the above processing and let B' be the 
corresponding beam set. Let dis{d') and B{d') be the original distance value and beam set at d'. 
Again, by the third algorithm invariant, this is the the first time C is processed in the direction 
from d to d'; hence, if dis{d') ^ oo, then d' must be labeled by a sweep from outside C and the 
beams of B{d') must enter C. Depending on the value of dis{d'), we may need to update the label 
of d' in several cases. 

1. If dis{d') = oo, we set dis{d') = 5 and B{d') = B'. Note that <5 > 2i + 3. Hence, the global-key 
of d' is strictly larger than that of d*, which has distance value 2i -|- 1. We insert d' into H (not 
Hr). 

2. If dis{d') < 2i-t-l, then since B{d*) = 0, the global-key of d' is no larger than that of d* regardless 
of whether B{d') is empty or not. By the first algorithm invariant, d' has been correctly labeled. 
We do a “post-processing” to compute the v-v-map in C by using the beams of Br{d) and B{d'). 
Again, we will give a corridor-post-processing algorithm for this step later in Section HI 

3. The remaining case is when dis{d') ^ oo and dis{d') > 2i -|- 1. By the same argument as in the 
previous case where B{d*) ^ 0, this case cannot happen. 

The above describes the right-sweep procedure. The left-sweep is similar. 

This finishes our discussion in the case in which B[d*) is empty for the diagonal d* at the root 
of H. As in the previous case, all algorithm invariants hold. 

The algorithm finishes if all three heaps H, Hr, and Hr become empty. After that, for each 
diagonal d in a junction cell, dis{d) and B{d) have been correctly computed. During the algorithm 
some corridors have been labeled correctly while others are left for post-processing. 

Specifically, consider any corridor C and let di and d 2 be its two doors with their beam sets 
B{di) and B{d 2 ). If C is not left for a post-processing, then C has been processed either from di 
to d 2 or from d 2 to di and the v-v-map in C has been computed after the processing. Suppose the 
above processing is from di and d 2 . Then, C is processed using the beams of B{di) and B{d 2 ) is 
obtained after the processing. We will show in Section H] that our corridor-processing algorithm on 
C runs in 0{m -I- {hi — /12 -|- 1) log hi) time, where m is the number of vertices of C, hi = |B(di)|, 
and /i 2 = |B(d 2 )|. If C is left for a post-processing, i.e., to compute the v-v-map in C by using 
B{di) and B{d 2 ), we will show in Section |4] that our corridor-post-processing algorithm runs in 
0{m -I- hi log hi /i 2 log / 12 ) time. 

Lemma 1. Given VD{B), our algorithm computes the v-v-map on VD{H) in 0{n-\-h\ogh) time. 

Proof. Recall that the beams in our algorithm are generated by locally-outmost diagonals. We say 
that two beams are different if they do not have the same generator. We say a diagonal is generated 
by a corridor if the generator of the diagonal is in the corridor. 

We first prove a claim that the number of different beams at the diagonals in all junction cells 
is 0{h). To see this, a key observation is that since a corridor is a simple rectilinear polygon, it 
can only generate at most two new beams that can go out of the corridor in the entire algorithm. 
Specifically, consider a corridor C with doors di and d 2 . Suppose the algorithm processes C in 
the direction from di to d 2 . Then, as will be seen later in Section HI if some beams of B{di) can 
directly illuminate d 2 , then all beams of B{d 2 ) are from B{di) (although some beams may become 
narrowed) and there is no new beam generated by C; otherwise all beams of B{di) terminate inside 
C and B{d 2 ) will only have one beam coming out of the corridor through d 2 (this is because C is a 
simple rectilinear polygon). Since any corridor can be processed at most twice, it can generate at 
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most two new beams. Since there are 0{h) corridors, the total number of beams on the diagonals 
of junction cells is 0{h). 

To obtain the running time of the algorithm, we analyze the time we spent on junction cells 
and the corridors separately. Due to the above claim, the size of B{d) for each diagonal d in any 
junction cell is 0{h). Since there are 0{h) diagonals in all junction cells, the total time we spent 
on processing them is 0{hlogh). 

For the time we spent on all corridors, it is the sum of the time of our corridor-processing 
algorithm and corridor-post-processing algorithm on all corridors. 

Let S denote the set of the diagonals in all junction cells. Define C, di, d 2 , B{di), B{d 2 ), m, hi, 
and /i 2 the same as before. Suppose C has been post-processed. Then the corridor-post-processing 
algorithm on C takes 0(m -|- hi log hi + /i 2 log / 12 ) time. The sum of the term m overall all corridors 
is 0{n). All beams of B{di) and B{d 2 ) are terminated inside C after the post-processing. Since the 
number of beams that are terminated inside corridors is no more than the number of different beams 
at the diagonals in all junction cells, by the above claim, the number of beams on the diagonals of 
S that are terminated inside corridors is 0{h). Hence, the sum of hi (resp., /i 2 ) over all corridors 
that have been post-processed is 0{h), and the sum of /iilog/ii -|- /i 2 log/i 2 over all corridors that 
have been post-processed is 0(/ilog h). Note that each corridor can be post-processed at most once. 
This proves that the total time of the corridor-post-processing algorithm in the entire algorithm is 
0{n + hlog h). 

We can use the similar approach to analyze the total time of the corridor-processing algorithm. 
Suppose C has been processed by the corridor-processing algorithm in the direction from di to d 2 - 
Then, the running time of the algorithm on C is 0{m + {hi — /12 + 1) log hi) time. Similarly, the 
sum of m over all corridors is 0(n). After the processing of C, the number of beams of B{di) that 
have been terminated in C is at least hi — h 2 . Since the the number of beams on the diagonals of 
S that are terminated inside corridors is 0{h), the sum of the term hi — /12 over all corridors that 
have been processed is 0{h). The sum of the additional term 1 over all corridors is clearly 0(h). 
Therefore, although a corridor may be processed twice, the total time of the corridor-processing 
algorithm in the entire algorithm is 0{n + hlogh). 

The lemma thus follows. □ 

The above computes the v-v-map on VD{B). Again, the above algorithm only labels diagonals. 
Using the similar approaches as discussed in Section 12.21 we can also label cells and maintain path 
information within the same running time asymptotically. We omit the details. 

The other three maps can be computed similarly. For computing the h-v-map on VD{B), one 
difference is on the initial steps, as follows. Initially, we let s generate two beams that are two 
horizontal rays towards right and left, respectively. We set the distance value of dg to 0, where dg 
is the vertical diagonal through s. Then, we consider dg as two duplicate diagonals associated with 
the above two beams respectively, and insert the two duplicate diagonals into the heap H. The 
remaining algorithm is the same as before except that we replace the distance values 2i -|- 1 and 
2i -|- 3 in the algorithm description with 2i and 2i + 2, respectively. The h-h-map and v-h-map on 
HD{B) can be computed in symmetric ways. 

We thus obtain the main result of this paper. 

Theorem 1. Given a set of h pairwise-disjoint rectilinear polygonal obstacles with a total of n 
vertices in the plane, after the free space is triangulated, we can construct a link distance map with 
respect to any given source point s in 0{n -\- hlog h) time, such that for any query point t, the link 
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Fig. 8. Illustrating VD{C): the two doors are shown with thick segments. 

distance of a min-link rectilinear s-t path vr can be computed in O(logn) time and vr can he output 
in additional time linear in its link distance. 

4 The Algorithms for Processing Corridors 

In this section, we present our corridor-processing algorithm and corridor-post-processing algorithm. 
In fact, as shown later, the second algorithm is simply to run the first one twice from the two doors. 
Below, we discuss the corridor-processing algorithm first. 

Let C be any corridor with two doors di and d 2 and their beam sets B{di) and B{d 2 ). Let m 
be the number of vertices of C, hi = \B{di)\, and /i 2 = \B{d 2 )\. Suppose we want to process C to 
compute the labels of all diagonals, by using the beams in B{di) and the distance value dis{di) 
of di- In particular, we want to obtain the beam set B{d 2 ) and the distance value dis{d 2 ) for d 2 . 
Our algorithm is conceptually similar to the DN algorithm if we apply it on C, but our approach is 
more efficient due to a better implementation by making use of the simplicity of the corridor (one 
crucial property is that there will be no merge operations on the beam sets). 

Denote by VD{C) the vertical decomposition of C (e.g., see Fig. [8]). Note that although di and 
^2 are now on the boundary of C, we still consider them as diagonals. Let C{di) denote the cell of 
VD{C) that contains di. Starting from C{di), we propagate the beams from B{di) to all other cells 
one by one until all diagonals have been labeled. 

Consider any diagonal d of VD{C) such that d is not di or d 2 . Since C is simply connected, d 
divides C into two polygons, and we use C{d) to refer to the one that does not contain di. Clearly, 
for any point t £ C{d), any path from di to t must intersect d. Consider any cell C of VD{C). The 
cell C may have diagonals on both its two vertical sides. It is not difficult to see that C has one 
and only one diagonal d such that C is in C{d), and we call that diagonal the entrance diagonal of 
C. Other diagonals of C are called exit diagonals of C. For example, in Fig. [8l d is the entrance 
diagonal of C, and d and d' are exit diagonals. Note that every diagonal is an entrance diagonal of 
one and only one cell. In particular, we consider di as the entrance diagonal of C{di). 

A general step of our algorithm works as follows. Consider a cell C and suppose the entrance 
diagonal d of C has been labeled, i.e., dis{d) and B{d) are available (and the beams of B{d) are 
stored in a balanced binary search tree T{d)). Initially, C is the cell C{di) and d is di. Our goal 
for processing d is to label all exit diagonals in C. 

Denote by ei and Cr the left and right sides of C, respectively. Without loss of generality, we 
assume d is on e/ (e.g., see Fig. [8]). Thus, the beams of B{d) are rightward. The distance values of 
the exit diagonals are easy to compute. If B{d) is empty, then we need to generate a single beam 
from the entire d and every exit diagonal of C obtains the distance value dis{d) + 2; otherwise, 
every exit diagonal of C obtains the distance value dis{d). Below, we focus on computing the beam 
sets of the exit diagonals. 
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Suppose ei has another diagonal d. Then, regardless of whether B{d) is empty, we set B{d) = 0 
(and T[d) = 0), which can be done in only constant time. 

Now consider any diagonal d' on e,.. First, we want to compute B[d) H d', i.e., the portions of 
the beams of B{d) that can illuminate d', and set B{d') = B{d) PI d'. If B{d) = 0, then since we 
generate a single beam from d, B{d') = B{d) r\ d' has at most one beam and can be computed in 
constant time. 

In the sequel we consider the case where B{d) ^ 0. Our main effort is on handling this case. 
Recall that the beams of B{d) are stored in a balanced binary search tree T{d). We add two special 
pointers to T{d) that point to the lowest beam and the highest beam of B{d) respectively so that 
we can access these beams in constant time. Note that with these special pointers, we can still 
perform the previous operations on T{d) each in logarithmic time. 

Depending on whether e,- has one or two diagonals, there are two cases. 

If er has only one diagonal d', let e be the obstacle edge on such that d' is the vertical 
extension of e (e.g., see Fig. ED. Clearly, er is the union of d' and e. Without loss of generality, 
assume e is lower than d'. 

Suppose 6 is a beam in B[d) and h' is the rightward projection of 6 on e,.. For any line segment 
I on Cr, we say that b intersects I if h' intersects I properly, and b fully intersects I if b' is contained 
in 1. 

Below we let b denote the lowest beam of B{d), which can be obtained in constant time by using 
the special pointers on T{d). We first check whether b intersects e. Depending on how b intersects 
e, there are three cases. The correctness of our setting in all these cases is based on that consists 
of d' and e from top to bottom (e.g., see Fig. ED- 

1. If 5 does not intersect e, then every beam of B{d) fully intersects d'. Thus we have B{d') = B{d) 
and T{d') = T{d). This can be done in constant time. 

2. If h intersects e but does not fully intersects e, then we set B{d') = B{d) and T{d') = T{d), but 
we also change 6’s length to that of its portion intersecting d'. This can be done in constant 
time. 

3. If b fully intersects e (so b “terminates” at e), then depending on whether B{d) has only one 
beam, there are further two subcases. 

If B{d) has only one beam, which is 5, then we simply set B{d') = 0 and T{d') = 0 since b 
terminates at e. 

Otherwise, we have \B{d)\ > 2. In this subcase, as before we use a split operation on T{d) to 
obtain T{d') and B{d') = B{d) D d' in 0(log |i?((i)|) time. Note that since b terminates at e, b 
is not in B{d'), and thus \B{d')\ < \B{d)\ — 1. This will be useful to our time analysis on the 
split operations. 

If Cr has another diagonal d", without loss of generality, we assume d' is the lower one. Let e be 
the obstacle edge on Cr- Then Cr consists of d", e, and d' from top to bottom (e.g., see Fig. fTop . 

Again, let h denote the lowest beam of B{d). Depending on how b intersects d', there are three 
cases. In every case below, we will also obtain a tree T' of beams that will be used later to compute 
the beams of d". 

1. If b does not intersect d', then we have B{d') = 0 and T{d) = 0. We set T' = T{d). 

2. If b intersects d' but does not fully intersects d', then B{d') consists of a single beam that is the 
portion of h intersecting d'. It is straightforward to construct T{d'). We set T' = T{d) but also 
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Fig. 9. The right side Cr of C consists of d! and e from Fig. 10. The right side Cr of C consists of d", e, and d' 
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change the length of b to that of the portion of b not intersecting d'. All above can be done in 
constant time. 

3. If b fully intersects d', then we further check whether the highest beam b* of B{d) intersects d'. 
If b* fully intersects d\ then we have B{d') = B{d) and T{d') = T{d). Also, T' = 0. 

If b* intersects d' but does not fully intersect d', then we set B{d) = B{d) and T{d') = T{d) 
but also change 6 *’s length to that of its portion intersecting d'. Also, we let T' only include 
the portion of b* not intersecting d'. All this can be done in constant time. 

If b* does not intersect d', then as before we use a split operation that split T{d) into two trees: 
T{d'), which consists of the beams of B{d') = B{d) D d', and T', which consists of the rest of 
the beams of B{d). This can be done in 0(log |i?((i)|) time. Let B' be the set of beams in T'. 
Note that since b fully intersects d', b is not in B', and since b* does not intersect d, b* is not 
in B{d'). Hence, we have \B'\ < \B{d) \ — 1 and \B{d')\ < |H((i)| — 1, which will be useful to our 
time analysis on the split operations. Further, since b fully intersects d' but b* does not intersect 
d', we have b' 7 ^ b* in this case, implying that |H((i)| > 2 . 

Next, we compute B{d'') and T{d") from the beams in the tree T', in the same way as the 
previous case where only contains one diagonal. Namely, we first check how the lowest beam of 
T' intersects e and then proceed accordingly for the three cases. 

The above labels all exit diagonals of the cell C. Based on our above algorithm, an easy obser¬ 
vation is that if \B{d)\ = 1, then \B{d')\ < 1 for each exit diagonal d' of C. 

Next, from each exit diagonal of C, we proceed with the same procedure. The algorithm is 
done once all diagonals have been labeled. The following lemma analyzes the running time of the 
algorithm. 

Lemma 2. Our corridor-processing algorithm on C runs in 0{m + {hi — /i 2 + 1) log hi) time. 

Proof. We make a few observations on our algorithm. First, for any diagonal d, if computing B{d) 
does not need a split operation, then it takes only constant time to do so. Second, for any diagonal 
d of C, \B{d)\ < |H((ii)| = hi always hold, and if \B{d)\ > 2, then the beams of B{d) are from 
B{di). Third, for any diagonal d, if B{d) is obtained by a split operation on a beam set B, then 
\B\ > 2 and \B{d)\ < \B\ - 1. 

The first two observations imply that if there are k split operations in the entire algorithm, then 
the total running time is 0{m + /cloghi). We claim that k = 0(hi — h 2 + 1). 

Indeed, by the second and the third observation, if a split operation is performed on B{d) 
for some diagonal d, then all beams of B{d) must be originally from B{di), and further, the split 
operation splits B{di) into two sets such that the number of beams in each set is at most |H(d)| — 1. 
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Therefore, it can be verified that if /12 > 2, then k = 0{hi — / 12 ); and otherwise, k = 0{hi). In 
either case, k = 0{hi — /i 2 + 1). 

We conclude that the algorithm runs in 0{m + {hi — /i 2 + 1) log hi) time. □ 

Next, we present our corridor-post-processing algorithm. Again, consider the corridor C as above, 
but now B{d 2 ) and dis{d 2 ) are also given as input and the beams of B{d 2 ) are towards the inside of 
C. Our goal is to compute the distance values for all diagonals of C by using the beams of B{di) and 
B{d 2 )- An easy solution is to use our above corridor-processing algorithm to process C twice, once 
only using dis{di) and B{di) and once only using dis{d 2 ) and B{d 2 )- Then, each diagonal d has 
been labeled twice, and we finally set dis{d) to the smaller distance value labeled above. Clearly, 
the total running time is bounded by 0{m + hi log hi + /i 2 log /i 2 ). 

5 Conclusions 

We presented a new algorithm for computing minimum link paths in rectilinear polygonal domains. 
The algorithm matches the f7(n -|- hlogh) time lower bound [23] if the triangulation of the free 
space can be done optimally in 0{n +hlogh) time. Our algorithm can also build the minimum-link 
distance map for a fixed source point to answer minimum link path queries. 

It would be interesting to see if our approach can be extended to work in a C-oriented world 
iH18|26j; one main stumbling block is that the DN algorithm does not generalize to C > 2 orien¬ 
tations. Another intriguing open problem is whether a result similar to ours, could be possible for 
paths and domains with unrestricted orientations m, e.g., whether there exists an 0{g{n) + f{h))- 
time algorithm for the problem where g is subquadratic (the 3SUM-hardness reduction [26| for the 
problem shows that / essentially cannot be subquadratic). 
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